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Abstract 

Current muon anomalous magnetic moment data have challenged Standard Model (SM) 
and seem to open a window for new physics. Since the difference between SM and experimental 
predictions is approximately 2.6(7. In the framework of the General Two Higgs Doublet Model 
(2HDM), we calculate the muon anomalous magnetic moment to get lower and upper bounds for 
the Flavour Changing (FC) Yukawa couplings in the leptonic sector. We also obtain lower bounds 
for the mass of the Pseudoscalar Higgs (m^o ) as a function of the parameters of the model. 
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Current muon anomalous magnetic moment data have challenged Standard Model 
(SM) and seem to open a window for new physics. Due to the 

high precision in value, it gives very restrictive bounds on physics beyond the SM. 
Although tte measurement is about 350 more precise is much more sensitive to New 

Physics since contributions to a; are usually proportional to mf. 

The most accurate measurement of hitherto, has been provided by the Brookhaven 
Alternating Gradient Syncrotron p|. Their data have an error one third that of the combined 
previous data f^, ref |0 reports 

= 11659202 (14) (6) x 10'^°. (1) 

On the other hand, SM predictions for has been estimated taking into account the 
contributions from QED, Hadronic loops and electroweak corrections. The final current 
result is [§ § 

af^ = 11659159.7 (6.7) x 10"^°. (2) 
Taking into account (0) is obtained 

Aa""/ = a;fP - af ^ = 42.6 (16.5) x 10"^°, (3) 

where a^^P is the world average experimental value |]T|. Consequently at 90% C.L. 

21.5 X 10"^° < AaJ"^ < 63.7 x 10^^°. (4) 

Aa^^ gives the room available for New Physics, so a^^^ differs from a^*^ approximately 
in 2.6cr. Therefore, physics beyond the SM is needed to achieve an acceptable theoretical 
experimental agreement. The most studied contributions to has been carried out in the 
framework of radiative muon mass models as well as the Minimal Supersymmetric Standard 
Model (MSSM), Eg string-inspired models, and extensions of MSSM with an extra singlet 
i- 

Moreover, a very interesting suggestion to conciliate the new experimental data with 
theoretical predictions is to consider models that includes FCNC at tree level. Interactions 
involving FCNC are forbidden at tree level in the SM, but could be present at one loop level 
as in the case of 6 — 57 0, — > ^~ 0, — K^ 0, t ^ 07 0] etc. Many extensions of 
the SM permit FCNC at tree level. For example, the introduction of new representations of 
fermions different from doublets produce them by means of the Z-coupling . Additionally, 
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they are generated at tree level by adding a second doublet to the SM [|^, such couplings 



can be gotten as well in SUSY theories without R-parity. Some other important new sources 
for FCNC might be provided by a muon collider, as the processes fi^ firmer) mediated 
by Higgs exchange [|rT|, |T2|, which produce Lepton Flavor Violation (LEV). 

However, there are several mechanisms to avoid FCNC at tree level. Glashow and Wein- 
berg |TB[ proposed a discrete symmetry to supress them in the Two Higgs Doublet Model 
(2HDM) which is the simplest one that exhibits these rare processes at tree level. There 
are two kinds of models which are phenomenologically plausible with the discrete symme- 
try imposed. In the model type I, one Higgs Doublet provides masses to the up-type and 
down-type quarks, simultaneously. In the model type II, one Higgs doublet gives masses to 
the up-type quarks and the other one to the down-type quarks. But the discrete symmetry 
T3| is not compulsory and both doublets may generate the masses of the quarks of up-type 



and down- type simultaneously, in such case we are in the model type HI [O]. It has been 



used to search for FCNC at tree level [jT5|, |[T6|| . 



Recently, the 2HDM type HI has been discussed and classified JT^, depending on how 
the basis for the vacuum expectation values (VEV) are chosen and according to the way in 
which the flavor mixing matrices are rotated. In brief, the reference [|17| shows that there 
are two types of rotations which generate four different lagrangians in the quark sector and 
two different ones in the leptonic sector. The well known 2HDM types I and II, could be 
generated from them in the limit in which the FC vertices vanish. It has been pointed out 
that the phenomenology of the 2HDM type HI is highly sensitive to the rotation used for 
the mixing matrices. 

In this paper, we calculate the contributions to Aa^^ coming from the 2HDM, which 
includes FCNC at tree level. We will constrain the FC vertex involving the second and third 
charged leptonic sector by using the result for Aa^^, equation (^). Additionally, we get 
lower bounds on the Pseudoscalar Higgs mass by taking into account the lower experimental 
value of Aa^^ at 90% CL and making reasonable assumptions on the FC vertex. 

The Yukawa's Lagrangian for the 2HDM type HI, is as follow 

-£y = 4Q,l^iU,r + vm^^iDjR + r]f~kL^iE,R 



+ ^^Q^L^2U,n + i^Q,L^2D,R + iffkL^2E,R + h.C. (5) 



where $1^2 are the Higgs doublets, r]ij and are non-diagonal 3x3 matrices and i, j are 



family indices. In this work, we are interested only in neutral currents in the leptonic sector. 
We consider a CP-conserving model in which both Higgs doublets acquire a VEV, 



(<^'i)c 





The neutral mass eigenstates are given by [|I8| 
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where tan/? = V2/V1 and a is the mixing angle of the CP-even neutral Higgs sector. Gz is 
the would-be Goldstone boson of Z and is the CP-odd neutral Higgs. 

Now, to convert the Lagrangian (§) into mass eigenstates we make the unitary transfor- 
mations 



El,r = {Vl,r) EIj^ 



from which we obtain the mass matrix 



m; 



diag 



Vj, 



V, 



R ; 



(8) 



(9) 



where M^"^ is the diagonal mass matrix for the three lepton families. From (|^) we can solve 
for obtaining 

^E,o ^ ^vlM^'^'Vn - ^77^'° (10) 

which we call a rotation of type I. Replacing it into (^, the expanded Lagrangian for the 
neutral leptonic sector is 

= -EM^^^'E ( sin aH^ + cos ah'') 

2Mwsmf3 ^ V ; 

2Mw 2Mw 

E7]^E [sin (a - (3) + cos (a - (3) 



n{I) 
'^Y(E) 



\/2sin/3 



Eri'^-f^EA" + h.c. 



(11) 



V2sm(3 

where the superindex (J) refers to the rotation type I. It is easy to check that Lagrangian 
(p!T|) is just the one in the 2HDM type I [|1^], plus some FC interactions. Therefore, we 



obtain the Lagrangian of the 2HDM type I from eq (|Tl|) by setting t]^ = 0. In this case it is 
clear that when tan /5 — > then rj^ should go to zero, in order to have a finite contribution 
for FCNC at tree level. 

On the other hand, from (Of) we can also solve for 7]^'^ instead of to get 



which we call a rotation of type II. Replacing it into the expanded Lagrangian for the 
neutral leptonic sector is 

- £^'p = -EM^'^'E ( cos aH'' - sin«/i°) 

^^^> 2Mwcos/3 ^ V J 

2Mw ^ 2Mw ^ 

+ E^^E [sin (a - p) + cos (a - P) 

v2cos/5 

+ J M '^l^EA' + h.c. (13) 
V 2 cos p 



The Lagrangian (0) coincides with the one of the 2IIDM type II ||T8[, plus some FC interac- 
tions. So, the Lagrangian of the 2HDM type II is obtained setting = 0. In this case it is 
clear that when tan (3 oo then E,^ should go to zero, in order to have a finite contribution 
for FCNC at tree level. 

In the present report, we calculate Aa^'^ in the 2IIDM with FC interactions. If we neglect 
the muon mass, the contribution at one loop from all Higgses is given by 



where 



^a^/ = ^Y.Fi^H.,m)ah (14) 



3 + m\ fmi - 4) + 2 In mi 
F{mH^.m) = ^3 (15) 



with rfiH- = mi/mij. and m; is the mass of the lepton running into the loop. The sum is 
over the index i = m/jo , m/^o , m^o . The coefficients are the Feynman rules for the FC 
couplings involved. 

If we take into account the experimental data (^), we get some lower and upper bounds 
on the mixing vertex ri{C,)^^ for the rotations of type I (II). In figure 1, we display lower 
and upper bounds for the FC vertices as a function of tan f3 for both types of rotations with 
rriho = tuho = 150 GeV and m^o oo. In the first case, rotation type I, the allowed region 



for rj^r is 0.07 < rj^-j- < 0.13 for large values of tan/3. Meanwhile, for rotation type II, the 
allowed region for small tan (3 is the same. From Lagrangian (11), which describes rotation 
type I, we can see that when tan/? 0, r^^^ should go to zero as well to mantain a finite 
contribution to Aa^. This behaviour can be seen from figure 1. For rotation type II occurs 
the same but in the limit tan (3 ^ oo. 

In figure 2, we show lower and upper bounds for the FC vertex as a function of m/^o for 
rotation of type II when m/jo = mno and m^o —>■ oo. We see that the larger value for tan/? 
the smaller value of . We only consider the case of rotation type II because there is a 
complementary behaviour between both rotations as could be seen in figure 1. In particular, 
for tan/3 = 1, the behaviour of the bounds for both rotations is the same. 

Observe that according to the Feynman rules from (0) and (p!3D, the Scalar (Pseu- 
doscalar) contribution to Aa^^ eq. (0,0) is positive (negative). Such fact permits us to 
impose lower bounds on the Pseudoscalar Higgs mass, by using the lower limit in eq. (^. 
According to this equation the room for new physics from g — 2 muon factor is positive 
definite, and it is a new feature from most updated results p[. 

Now, to take into account the experimental value (^), we should make a supposition 
about the value of the FC vertex. A reasonable assumption consists of taking the geometric 
average of the Yukawa couplings [0] i.e. t] (O^r ~ 2. 5 x lO^'^. Additionally, we shall use also 
the values t] (O^r ~ 2. 5 x 10~^ and rj (O^r ~ 2. 5 x 10~^ which are one order of magnitude 
larger and smaller than the former. Using these suppositions and the experimental value 
(I) we get lower bounds for m^o and they are plotted in figures (3)- (5). 

Figure 3 displays rriAO vs tan/3 using rotation type II with the three values of 
C,^r mentioned above and setting m/^o = uiho with m/jO = 110, 300 GeV. It could be seen 
that in the limit of large tan/3, the lower limit reduces to rriAO ~ rnf^o. The same behavior 
can be seen in rotation type I but the bound m^o ^ mi^o is gotten in the limit of small tan /3. 
We also see that the smaller value of the stronger lower limit for tuao ■ 

Figure 4 shows m^o vs m/^o with = 2.5 x 10"'^, 2.5 x 10^^ and tan/3 = 1, using 
m^o = rriHO and m^o = 110 GeV. With this settings, the value = 2.5 x 10~^ is excluded. 
Using such specific arrangements, the bounds are identical in both types of rotations. 

In figure 5 we suppose that mf^o = 110 GeV, 171^0 = 300 GeV. The figure above shows 
the sensitivity of lower bounds on m^o with the mixing angle a, for rotation type II, taking 
tan/3 = 1. The value = 2.5 x 10"'^ is excluded again. The constraints are very sensitive 
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to the a mixing angle for t,^n = 2.5 x 10~^ but rather insensitive for ^^t- = 2.5 x 10~^. The 
figure below shows m^o vs tan/3 for m^o = 110 GeV, ttlho = 300 GeV, a = n/G, for rotation 
type II and considering the same three values of C,^r- The m^o lower asymptotic limit for 
large tan/5 is approximately rriho. 

In conclusion, we have found lower and upper bounds for the FC vertex t] (O^r 
text of the general 2HDM by using the allowed range for Aa^'^ at 90% CL and utilizing sev- 
eral sets of values for the parameters of the model. Additionally, in the limit m^o —>■ oo, we 
get that for small (large) values of tan /? the allowed range for the FC vertex rj^^ (^^r) becomes 
narrower, and both upper and lower bounds go to zero in the rotation of type I (II). 

On the other hand, we have gotten lower bounds on the Pseudoscalar Higgs mass of the 
2HDM coming from the g — 2 muon factor, by using the experimental value of Aa^^ and 
making reasonable assumptions on the FC vertex ?7(0^r- Specifically, we have taken for 
r] iO/iT geometric average of the Yukawa couplings, and we also utilized values one order 
of magnitude larger and one order of magnitude smaller. Taking these three values for the FC 
vertex we find that the smaller value for rj {^)^^ the more stringent lower bounds for m^o . 
Additionally, assuming rriHO = mhO, we show that in the limit of small (large) tan/3 the 
lower bound of m^o becomes merely m^o ~ rriho for rotation of type I (II). In the case of 
different scalar masses, there is still a lower asymptotic limit for m^o. Notwithstanding, 
these lower constraints on rriAO should be consider carefully, since for rj (O^r '^^'^ only 
make reasonable estimations but they are unknown so far. 
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FIG. 1: Lower and upper bounds for rj^r (Cmt) vs tan/?, for rotations I and II using m/jO = m^o 
150 GeV and m^o — oo. 
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FIG. 2: Lower and upper bounds for vs m^^o, for rotation of type II, taking m^o = m^o and 
rriAo —>■ oo, the pair of short dashed hnes correspond to tan/3 = 0.1, the long dashed hnes are for 
tan/3 = 1, and the sohd hnes are for tan/3 = 30. 
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FIG. 3: Contour plot of m^o vs tan/3 using rotation type II and assuming m/jO = rrif^o. Short 
dashed hues correspond to = 2.5 x 10~^ for m^o = 110 GeV (below) and m^o = 300 GeV 
(above). Long dashed lines correspond to = 2.5 x 10~^ for m^o = 110 GeV (below) and 
mjjo = 300 GeV (above). Finally, solid lines correspond to = 2.5 x 10~^ for nifjo = 110 GeV 
(below) and mfjo = 300 GeV (above). 
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FIG. 4: Contour plot of m^o vs m^o setting tan/3 = 1, Long dashed lines correspond to 
2.5 X 10~^ for m/jO = 110 GeV (rigth) and m^o = m^o (left). Solid lines correspond to 
2.5 X 10"^ for m/jO = 110 GeV (below) and m/jO = nifjo (above). 



12 



1000 



800 



600 



40Q 



200 




lOQO 
800 
600 
400 
200 


5 10 15 20 25 30 

FIG. 5: (top) Contour plot of m^o vs a, for rotation type II, with m^jO = 110 GeV, nifjo = 300 
GeV and tan/3 = 1. Dashed hne correspond to = 2.5 x 10~^ , solid line correspond to 
^fj,r = 2.5 X 10~^. (bottom) Contour plot of m^o vs tan/3 for ruf^o = 110 GeV, m^^o = 300 GeV, 
a = 7r/6, and for rotation type II. Short dashed line correspond to = 2.5 x 10~^, long dashed 
line correspond to = 2.5 x 10~^, and solid line correspond to = 2.5 x 10~^. 
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